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Abstract 
The combined effects of climate and density-dependence on ungulate population were 
widely detected, which played the critical role in maintaining sustainability of 
rangeland ecosystems. But, previous studies in Mongolia plateau seldom considered 
the density dependent effect of the domestic animals, but mainly focused on the 
impacts of natural hazards. In order to fill this knowledge gap, we selected ten 
long-term sheep population series from the eastern part of Mongolia plateau and 
detected the density dependent effects. Three kinds of methods were independently 
employed in this study, i.e. the demographic parameters analysis, statistic tests, and 
the model selection method. 

The density dependence was detected by all methods during the nomadic period 
of 1962-1979. It found that three populations were density dependence. The seven 
other populations were only regulated by density dependent effect during normal 


winters, and just controlled by climate during harsh winters. The limitation of 


density dependent effect can explain the reason of why herdsmen migration distance 
positively correlated with the size of their herds. “Therefore, there were two types of 
natural pushing force to drive nomadic people migration, one was the natural hazard 
which was widely accepted, and the other was the density dependent effect. 
Especially, “pushing off of density dependence” can explain the reason why Genghis 
Khan’ rise and conquest happened during the medieval climatic optimum, when the 
natural hazard was less than other period in the record history of Eurasian. 

On the other hand, the study showed that the climate effect on population had 
weakened and the density dependent effects disappeared during the sedentary period. 
This finding is the reason of why the sheep population had showed a continued 
growth trend and caused “overgrazing” problem. Our results indicated that the 
disappearing of natural negative feedbacks had driven the ecosystem regime to shift 
from the “animal population is low but ecosystem is spontaneously sustainable status” 
into “animal increased rapidly but vegetation degenerated status” during the social 
change of sedentary in the Mongolian plateau. 
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Introduction 

Many nomadic empires had suddenly risen in Eurasian Steppe. “They attacked the 
European and Asian sedentary countries frequently; transmitted the languages, 
techniques, religions and plagues between eastern and western Eurasia through the 
Eurasian steppe corridor (Chase-Dunn et al. 2010). Especially, the rise and conquest 
of Genghis Khan and his empire were the most outstanding events in Eurasian steppe 
history which had shaped the world (Turnbull 2003). Many scholars believed that 
climate change may be a reason of the rise of nomadic empires, and have developed 
many hypotheses to explain it. Earlier scholars considered the drought and climate 
cooling may cause the nomadic invasions by destroying the domestic animal 
population and bringing subsequent famine (Huntington 1907, Toynbee 1934, Jenkins 
1974, Lamb 1997, Brown 2001). But, recent climate reconstruction data showed the 
climate in Genghis Khan’s rise and heyday were warmer and wetter than normal the 
past 900 years (HVIS 2011), and there is no evidence of climatic sudden cooling 
before the Mongolian conquest (Wu et al. 2009). Scholars recently suggested that 
the vegetation and domestic animal should flourish during the Medieval Climate 
Optimum, which assisted the rise and conquest of Mongolian (Wt et al. 2009, HVIS 
2011). However, this climate optimum hypothesis can not explain what forced the 
nomadic people to migrate out if the climate disasters were absent from natural 
perspective. Compared with the abundant quantitative studies about climate 
reconstruction data and historical records, the logic link between climate change and 


human migration was usually accepted as the common sense rather than analyzed by 


science evidence (Perch-Nielsen et al. 2008). The common sense was that only 
climate disasters can drive nomadic people to migrate out steppe though the collapse 
of nomadic animal population. Hence, in order to fix the weak link of these 
hypotheses, this paper would examine other nomadic animal population regulation 
mechanisms in Eurasian steppe. 

Mongolia plateau is the last nomadic region in the world (Sneath 1998, 
Humphrey and Sneath 1999), and also a representative part of Eurasian Steppe (Fig. 
1). Thus, the knowledge from the nomadic animal population dynamics in nomadic 
Mongolian grassland can help us to understand the nomadic history of Eurasian 
steppe. Both wild and nomadic domestic ungulates populations usually experience a 
“boom and bust” cycle in Mongolia plateau (Begzsuren et al. 2004, 
Fernandez-Gimenez et al. 2012), and the other part of Eurasian Steppe, such as in 
Kazakhstan grassland (Robinson et al. 2003). The ungulate population crash was 
usually caused by many kinds of climate hazard. Especially, blizzard and following 
snow-locked pasture in winter were the most serious disasters in Eurasian steppes, 
which were named as “duzd” in Mongolian traditional ecological knowledge 
(Begzsuren et al. 2004, Fernandez-Gimenez et al. 2012). Impacts of dzud on 
ungulates population and human livelihood were so serious that previous studies 
usually suggested that climate disaster dominate the ungulates population dynamics in 
Mongolian plateau (Begzsuren et al. 2004, Vetter 2005, Zhang and Li 2008, Li and 
Huntsinger 2011, Li and Li 2012). 


However, population ecological studies found that animal population dynamics 


were influenced by the combined effects of intrinsic factors (e.g. density-dependence) 
and extrinsic factors (e.g. climate factors, predation, and harvest); especially 
density-dependence and climate factors can regulate ungulates population dynamics 
without predator and harvest in both mild and harsh environments (Sæther 1997, 
Forchhammer et al. 1998, Post and Stenseth 1999, Aanes et al. 2000, Forchhammer et 
al. 2001, Solberg et al. 2001, Jacobson et al. 2004, Lima and Berryman 2006, Tyler et 
al. 2008, Colchero et al. 2009, Simard et al. 2010, Herrando-Perez et al. 2012, Imperio 
et al. 2012). Some domestic ungulates population dynamics are also controlled by 
density dependence and climate factors just as their wild neighbors in cold arctic 
region or dry African grasslands (Desta and Coppock 2002, Weladji and Holand 2003, 
Helle and Kojola 2008), where the environment may be as harsh as Mongolia plateau. 
Now, only very few studies mentioned the impacts of density dependence on 
domestic animal population in Mongolia plateau. Neglecting density dependent 
effects on ungulate population dynamics may cause some knowledge gaps in 
understanding the rangeland ecosystem regulation mechanisms. At least, current 
knowledge can not tell us how the ecosystem regulates the animal population under 
the climate optimum condition in Eurasian steppe. 

On the other hand, understanding ecosystem regulation mechanisms on domestic 
animal population has practical value for current ecosystem managements. The 
herdsmen society in Inner Mongolia had transform from nomadic style into the 
sedentary one (Sheeny 1993, Sneath 1998, Humphrey and Sneath 1999), and this 


social transition is happening in Mongolia now (Johnson et al. 2006). After the 


social transition, the animal population showed a rapidly increased trend rather than 
the bloom and bust cycle pattern (Li et al. 2012), and then overgrazing issue had 
emerged which caused serious environment problems, such as vegetation denegation, 
desertification and biodiversity loss (Tong et al. 2004, Jiang et al. 2006, Huang et al. 
2009, Li et al. 2012). Because the rangeland ecosystems resilience was maintained 
by the negative feedback of population density and limitation of climate, which can 
prevent ungulates to overexploit vegetation (Forchhammer et al. 1998), so analysis 
the change of rangeland ecosystem feedback mechanism can explain the reason why 
the steppe resilience weakened after the social transition, and provide the advice for 
future rangeland ecosystem sustainable managements. 
In this study, we collected ten long-term sheep population time series data from 
a eastern part of Eurasian steppe to analysis the ecosystem regulation mechanism of 
ungulate population in nomadic and sedentary periods (Fig. 1). Sheep is the major 
domestic species in Eurasian steppe, composites nearly 60% of total domestic animal 
population in Mongolia plateau, and is the important economic base of nomadic 
society and herdsmen livelihood(Sternberg 2008, Zhang and Li 2008, Li and 
Huntsinger 2011, Li and Li 2012). Many influential ecologists have suggest the 
animal population growth can be separated into density dependent component, density 
independent component, and “the component that represents the losses from harvest” 
(Beddington and May 1977, Leigh 1981, May 1981 , Lande et al. 1995, Saether et al. 
1996, Engen et al. 1997, Kokko and Lindstrom 1998, Diserud 2001). This 


separation can be written as following (Saether et al. 1996): 


AN SPIN =F CV AN (1) 
Where AN is the annual increase population, r(t)N is the density-independent 
growth number, /(N)is the losses due to density dependent effect on population, 
and A(N) is the losses due to harvesting. Using this method, we can separate the 
harvesting effect from intrinsic growth rate of domestic population, and focus on 
density dependence and climate effect on domestic sheep population. This study 
would answer following scientific questions: 1st, whether the sheep population is 
regulated by the density dependent effect in Mongolia steppe? 2nd, what are the 
relationships between the impacts of climatic disaster and density dependence on 
sheep population dynamics? 3rd, can these findings of this study fix the weak link 
of current hypotheses about the nomadic empire rises and conquest? 4th, did the 
effects of density dependent and climate on sheep population dynamics changed for 
the social transition (from nomadic style to a sedentary ones), and what the influence 


of this social changes on ecosystem resilience? 


Method 


Study area 
East Ujimugin grassland is located in the eastern part of Eurasian steppe, covering a 
total area of 47,000 Km? and ranging from 44.7? to 46.6? in latitude and from 115.1? 


to119.7? in longitude (Fig.1). The winter is long and cold. During the whole study 


period (1961-2005), the mean annual temperature was -18.7°C in winter (January, 
February and December in the yesteryear), and -1.2 °C in early spring (March and 
April); the daily temperature can be below -40 °C in extreme cold days. In the study 
period, mean annual winter precipitation was 13mm, and the extreme value was 
40mm in 1978. But Imm precipitation can form 15mm depth snow, and the snow 
does not melt during the whole winter when the daily high temperature is below -7°C. 
Thus if the snow depth reach a threshold, the pasture is locked until the next spring 
coming. Both single blizzard and several snow events can form the snow disaster as 
long as snow depth reaches the threshold. But this threshold is site-specific and 
influenced by plant height and topography. The major vegetation types are typical 
steppe and meadow steppe in East Ujimuqin grassland. There are also some 
“non-zonal” vegetation types such as shrub lands and sandy lands. The primary 
productivity of the steppe ecosystems peaks in August and is determined by rainfall in 
growing season (Wu and Loucks 1992, Li et. al. 2012). There are many lakes, 
springs, and rivers in this grassland, thus the animal can drink water easily in the 
growing season. 

This region is run by a Mongolian autonomy government during the whole study 
period. Almost all the rural residents are Mongolian herdsmen, who keep intact 
Mongolian traditional culture and ecological knowledge. Ujimuqin sheep is a 
native species which is domesticated by Mongolian ancestors. Before 1980s, local 
herdsmen live a traditional nomadic life, sheep are exposed in natural environment 


all the year round (Li and Ma 1993, Wang 2006, Jia 2011). After 1980s, the 


Mongolian residents had gradually settle down (Li and Ma 1993, Sheeny 1993). 
During sedentary period, herdsmen let sheep live outdoor in sunny days of winter, 
but they started to store forage to feed their animal in snow days. The social 
transition from nomadic life and sedentary way was not suddenly happened. Thus, 
we define the three types of landuse period in this research, i.e. traditional nomadic 
period from 1961 to 1979, sedentary period from 1987 to 2005, and social transition 
period from 1980 to 1986. We mainly studied the sheep population dynamics in 


nomadic and sedentary period in this research. 


Data acquisition and pre-processing 

We chose ten continuous sheep census data in this study, and divided them into 
nomadic period and sedentary period. The annual sheep census data between 1961 
and 2005 was acquired from East Ujimuqin agro-pasture management bureau. The 
data sets contain annual population (N), birth number (B), mortality number (M) and 
harvest number (H, the sum of slaughtered and sold numbers). The data are 
censused at sumu level in every June. Sumu is composed by three or four villages, 
and cover an average area of 2000Km” (Fig 1). Mongolia herdsmen do not eat 
and sell dead animal, thus the mortality number are independent with harvest number. 
The herdsmen mainly slaughtered and sold their animal after annual census, and adult 
females animal are seldom been harvested. “Therefore, harvested individuals can not 
reproduce in next census year. Thus, the relation among each parameter can 


describe as following: 


Besides, we also chose an integrated climate parameter to describe the winter climate 
condition, which was named as winter index in following. “This winter index was 
calculated from the four critical winter climate parameters which were recorded by 
China meteorological administration. The four single climate parameters were the 
winter snow depth, the pasture lock-out days, the extreme cold days in cold season, 
and the average temperature of reproduction period (early spring). This integrated 
climate parameter can represent climate condition in multiple periods such as winter 
and reproductive period, and can describe average and extreme simultaneously, thus it 
is better than these single climate parameters in explaining the population dynamics. 
The calculation and discuss of the winter index were recorded in the support 


information]. 


Estimating impacts of climate and density dependence from demographic 
parameters 

At first, we described the sheep population and demographic parameters dynamics in 
nomadic and sedentary periods respectively. “Then, we investigated the effects of 
climate and population density on birthrate, mortality by regression analysis in the 
two landuse period at each site. “These regression analyses used the regular linear 
equations, regular non-linear equations and threshold equations to fit the models. 
And then we used the coefficient of determination (R”) as criterion to select the better 


regression relationship. In order to display regressions pithily, we plotted these 


regressions at the region scale. Before the plotting, we standardized each time series 


to remove site-specific noise. The standardization method is as following: 
WEAN j RENN TE AY PA WE FA PA ————————" 3) 


where the x, xs stand for the original data and standardization data respectively, xis 


the average of a time series, and s is the standard deviation. This method transform 
the original time series into standardization ones, which average was 0 and standard 


deviation was 1. 


Statistic tests for density dependence 

We employed three representatively statistic methods to detect the density dependence 
effect from the 10 sheep population time series in two landuse periods. The three 
methods are Bulmer R(R*) test (Bulmer 1975), Pollard et al.'s test (Pollard et al. 1987) 
and Dennis and Taper's test (Dennis and Taper 1994). All the null hypothesis of the 
three tests were density independence. But the statistics which were employed in 
these statistical inferences did not follow regular distribution. Thus, these tests used 
random shuffle method and bootstrap method to generate statistic distribution 
respectively. Besides, all three methods were adjusted to adapt the harvesting 
scenario in our studying. “The adjusting followed the equation 1, and based the fact 
that the harvest component does not contribute the population growth in the next year. 
Besides, Bulmer R (R*) test are less effective when the population data with trend 
(Shenk et al. 1998, Jacobson et al. 2004), hence we did not use this method to test the 


data in sedentary period which showed significant growth trend. The details of these 
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statistic tests are in supporting information 2. 


Model evaluation for density dependence 

We also used regression models and model selection method to estimate winter 
climate and population density effects on sheep population dynamics. “The harvested 
effects on population were separated out. Thus, the intrinsic growth rate is defined 


as following: 


MON eX E metn) "————— — (4) 


t 


We employed regular linear models and threshold models to describe effects of 
climate and population density on intrinsic growth rate, and used model selection 
method to find the optimal model for each site at the two landuse periods. “The 
regular regressions were started from a whole model which included the terms of 
winter climate, population density and the interaction term of them. Then, the others 
regular regressions were derived from the whole model by reducing different terms. 
We also induced many kinds of threshold and generated several threshold regression 
models. At last, the correction of Akaike information criterion (AICc) was 
employed to select the optimal model. The details of this part can found in 
supporting information 3. We also plotted the optimal regression at the region scale 


to display the relationship clearly. 


Simulating the population dynamics of 1987-2005 under the nomadic scenario 


We used the real population data in 1987 and the climate parameter in the period of 
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1987-2005 to simulate the population dynamic under the nomadic natural-society 
scenario. The natural population regulation mechanisms in the nomadic period were 
described by the optimal population model of nomadic period at each site. In order 


to describe nomadic scenario, the harvest rates of each site were defined as follows: 


Where H is the mean harvested rate of each site at nomadic period, S is standard 


deviation, and ris random number which in the interval of [-1, 1]. “Then we could 
compare the differences between real animal population demography and simulation 


ones in 1987-2005. 


Result 


Demographic patterns 

All sheep population fluctuated dramatically in nomadic period (1962-1979). The 
winter climate was mild in 1960s and harsh in 1970s (Fig. 1 of supporting information 
1). Sheep population had increased very fast in early period of 1960s, but the 
increase speed gradually reduced before the snow disaster winter of 1969-1970 when 
all ten populations crashed in this year (Fig. 1). After the population collapse of 
1970, all population recovered rapidly, then the population growth was slowdown 
again, at last the huge duzd winter of 1977-1978 caused another regional population 


crash (Fig. 1). However, in sedentary period (1987-2005), sheep population showed 


a continued increased trend in all the sites and doubled in most sites. Besides, sheep 
population seldom crashed during harsh winters of sedentary period, i.e. 1990-1993 
and 2001 (Fig. 1). 

In nomadic period, demographic parameters reached the peak of mortality and 
the nadir of birthrate during disaster years. Besides, the mortality showed increased 
trend and the birthrate decreased as the population increased in these no snow disaster 
years of this period (Fig 2). Impacts of duzd on mortality and birthrate still exist 
during sedentary periods, but the morality was constant and birthrate was increasing 
in the duzd-absent years of this periods. Generally speaking, the mortality was 
higher in nomadic period than sedentary period. At regional level, the average of 
mortality in nomadic period was 0.06 with the variance of 0.07, while the value was 
0.02+0.03 in sedentary period (Fig. 2). But, the birthrate was lower in nomadic 
period (0.41+0.12) than sedentary period in which the value is (0.56+0.14). At last, 
the human influence in nomadic period was weaker than sedentary period. The 
mean harvested rate was 0.22+0.10 in nomadic period. But this number was 
0.28+0.19 in the period of 1987-2000, and increased to 0.43+0.07 in the period of 


2001-2005 (Fig 2). 


Effects winter climate and population density on mortality and birthrate 
In the nomadic period, the mainly climatic reason for mortality fluctuation was winter 
snow condition changes, for example the snow depth can explain 53% of total 


mortality variance at regional level. But the birthrate was primary impacted by 


extreme cold events among all climatic variables, which can explain 29% of total 
birthrate variance at regional level in the same period. The winter index, which was 
integrated from four climate parameters, can explain 63% of total mortality variance 
and 32% of total birthrate variance in nomadic period. Thus, the winter index was 
more effective than single climate parameters and employed in following study. The 
details of results and discuss about climatic parameters selection were represented in 
Supporting information 1. 

The relationships between winter index and the demographic parameters show a 
threshold trend in nomadic period. Especially, the linear threshold relations are clear 
at the regional scale (Fig. 3). When the winter index were less than mean value plus 
1/2 the standard error of winter index (Mean + 1/2 SD), the R? of regressions between 
winter index and demographic parameters were very low at regional scale (for 
mortality, R7=0.00, p=0.7876, for birthrate and R7=0.06, p=0.5286, Fig. 3); otherwise, 
all regressions are significant and with high R? (for mortality, R7=0.77, p«0.0001; and 
for birthrate, R?=0.54, p«0.0001, Fig. 3). We defined this value as a climatic 
threshold for following analyses in this study; when the annual winter index was 
under the threshold, the winter was defined as mild; otherwise, the year was defined 
as harsh winter. 

In nomadic period, the relation between population density and demographic 
parameters still followed a threshold nonlinear shape. Regular linear regressions 
were no significant at most sites in the whole nomadic period (Table 1 of supporting 


information 4). However, it found that population density only significantly 


correlated with mortality and birthrate at the mild years. In harsh years, the 
regressions between population density and demographic parameters were not 
significant at each site (Table 3 of supporting information 4). At the regional scale, 
this piecewise phenomenon can be found easily, the regression population density 
against mortality and birthrate were both significant in mild winters, (for mortality, 
R?=0.22, p«0.0001; and for birthrate R?-0.24, p«0.0001; Fig. 4). But the two 
regressions were not significant in harsh year (for mortality, R7=0.22, p«0.0001; and 
for birthrate R220.24, p«0.0001; Fig. 4). 

However, the effects of climate and population density on demographic 
parameters weakened in sedentary period of 1987-2005. The winter index can 
explain 56% of total mortality variance and 14% of total birthrate variance at regional 
level in sedentary period, which are lower than nomadic periods. On the other hand, 
the regressions between population density and demographic parameters in most sites 
were not significant in this period. Besides we did not find any threshold evidences 


from these relationships in the sedentary (Supporting information 4). 


Density dependence evidence from Statistic results 

Bulmer R test, Pollard et al.’test and Dannis and Taper’ test suggested density 
dependence only exist in a few population. In nomadic periods, all the three statistic 
tests suggested the density dependence in the population of sitel and site7 (Table 1, 2 
and 3). The population in site6 was proved as density dependence by Pollard et.al’ 


test and Dannis and Taper test respectively (Table 2 and 3), besides the Bulmer R 


Statistic value of site 6 population was slightly larger than a=0.05 critical value. 

Thus, we considered the site 6 population was density dependence too. The statistics 
values, which indicated the feedback of population density in Pollard et. al’ test and 
Dannis and Taper’s test, were negative for population of site 1, 6 and 7. These 
statistics were named “Tax” and “b”, the calculation method of the two statistics was 
presented in supporting information 2. In sedentary period, Pollard et al.’ test and 
Dannis and Taper’s test suggested the population were density dependence of site 4 
and site 5 (Table 2 and 3). But all the rax and b were positive which reflected the 
positive feedback of population density for the population of the two sites in the 


period of 1987-2005. 


Model selection results 

We observed the following nonlinear phenomenon in nomadic periods: there were 
significant linear correlation between population density and intrinsic growth rate in 
mild years, but the linear correlations were not significant at harsh years for most sites 
(except sitel, 6 and 7). “We also presented the phenomenon at regional level (Fig5). 
The R? of regression between population density and intrinsic growth rate was 0.35 
(p<0.0001) in mild years and 0.03 (p=0.31) at harsh years. On the other hand, the 
climatic effect was very weak in mild years, but very strong in harsh years. “The 
regional R? of regression between winter index and intrinsic growth rate was 0.09 
(p=0.004) in mild years and 0.85 (p«0.0001) in harsh years (Fig 5). Thus, we 


presented a candidate threshold model: 


y; = 


a — cd, ME AAS < nee kn koko ke kk TRET epa REL Oke kaa HEE QE Era (6) 
b—dw; if — w; Zihreshold. 


This threshold model means the intrinsic growth rate was just controlled by 
population density in mild years, and just regulated by climate in harsh winter years. 
It also indicated that the effects of climate and population density on sheep population 
dynamics followed a linear threshold form but without interaction. “This type of 
threshold model had the smallest AICc and the highest R? among all regular and 
threshold models for all sites expect sitel and site6 (Table 4). On the other hand, the 
best models of sitel and site6 were regular linear models without the interaction term. 
All coefficients of climate term or population density term were negative in the 
optimal models in site 1 and site 6, which means the sheep population was density 
dependent in all winter condition (Table 4). At last, model efficacies of regular and 
threshold model about population of site 7 were not significant different with each 
other, for the difference of AICc was less than 1. So, the population of site 7 still 
can be considered as density dependence. 

In sedentary period, the model selection results suggested the negative effects of 
winter climate and population density on population growth had weakened or 
disappeared. All threshold models were excluded by the AICc criterion and their 
explanation efficacies were lower than regular models which meant the nonlinear 
phenomenon disappeared. “The optimal models for each site can be divided as three 
types: the intrinsic growth rate was just influenced by climate, just influenced by 
population density, or influenced by combined effects climate and population density 


without interaction (Table 5). But, all the R? of optimal models were relatively low. 


However, all the coefficients of population density term were positive, which means 


the population abundant had sped population growth in sedentary period. 


Simulated demographic pattern 

The simulations suggested that the demographic pattern were boom and bust in all 
sites under the nomadic social-ecosystem scenario with the real winter climate data of 
the period 1987-2005, which were different with the actual population dynamics (Fig. 
6). Besides, simulations population abundant data did not exceed their peak value of 
nomadic abundant in all sites. Therefore, the reason of why the population had 
continued growing in sedentary period was that the limitations of climate and density 


dependence had weakened or disappeared. 


Discussion 


Density effect did exist in nomadic period (1962-1979) 

This study proved that the density dependence was an important regulation factor on 
sheep population dynamics in nomadic period of Mongolia grassland. Tn our study, 
some sheep populations were density dependent, the other populations were also 
regulated by population density at non-duzd winters. But, only a previous study had 
the similar finding with this study in Mongolia plateau, which found that there were 


positive correlations between yesteryear population density and mortality during duzd 


year (Tachiiri et al. 2008). So we would still discuss our results cautiously in this 
section. 

At first, there are some weaknesses in the statistic methods; for example only 
Bulmer R* test is not sensitive against census error among the all statistic methods 
(Fox and Ridsdillsmith 1995), but efficacy of Bulmer R* was lowest among all 
statistic methods (Shenk et al. 1998). Freckleton et al. (2006) suggested that two 
common types of census error can reduce the statistic efficacies. “These census errors 
are usually caused by the two following reasons: the population abundant were not 
been fully censused, but just estimated by few individual; or the objective population 
was not a full self-maintain population rather than a subset population (Freckleton et 
al. 2006). But these error types did not exist in the census data of our study, which 
recorded all animals in a geographic region. The error of our census data maybe root 
in the political and economic reasons. Local officers often aggrandized the 
population increase abundant to gain political achievement during mild years, which 
underestimated the negative effect of population density; they would also overdrew 
the population loss during duzd year to get more alms from central government, which 
overestimated the impact of snow disaster. The political and economic census error 
increased the probability of making type II error of these statistic methods (accepting 
the null hypothesis of density independence). Besides, the model selection results 
confirmed these statistic results. “The optimal models were regular linear models 
with negative density term for these three populations of site 1, 6 and 7, which means 


population growth was limited by the population density. Hence, these findings of 
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density dependence were robust. 

On the other hand, for these density independence populations which were tested 
by statistic methods, population density still impacted the population dynamics during 
non-duzd years. “These optimal models of these density independent populations 
were linear threshold equations. Besides. the birthrate and mortality showed 
significant nonlinear threshold forms in responding climate and population density 
changes too. “These threshold evidences indicated that the population density was 
major driven factor of population dynamics at mild winter years, while only climate 
impacted population dynamics at harsh winter years. Besides, it was common 
knowledge in ungulate population ecological studies that the density independence 
were caused by the nonlinear or threshold influences of population density and 
climate (Coulson et al. 2001, Mysterud et al. 2001, Jacobson et al. 2004, Chan et al. 
2005, Lima and Berryman 2006, Tyler et al. 2008, Colchero et al. 2009). Because 
pervious studies mainly focused on the duzd event in Mongolia plateau (Begzsuren et 
al. 2004), thus these studies had very few chance to detect the impacts of density 
dependence. Therefore, both population density and climate can regulate the 
domestic animal population dynamics in nomadic period of Mongolia plateau. 

To our surprise, we did not detect the interaction effect of climate and density 
dependence on population dynamics; for all optimal models did not contain 
interaction term of population density and climate. Although numerous of studies 
found the interaction effect of climate and density dependence (Aanes et al. 2000, 


Forchhammer et al. 2001, Weladji and Holand 2003, Jacobson et al. 2004, Lima and 
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Berryman 2006, Tyler et al. 2008), while there were also some previous studies did 
not discover the interaction effects (Solberg et al. 2001). The interaction effects of 
climate and population density usually root in following reasons. The first reason 
was that the population density aggravates the food shortages by competition during 
harsh climate periods. “The second reason was that snow can cause positive effects 
on animal population growth by improving the quality and quantity of forage in the 
next growing season; many previous studies had found the effects in Europe 
(Mysterud et al. 2001). But, the food competition hardly happened during the 
extreme harsh winters. Because all plants were covered by snow, thus there was no 
food for competition. On the other hand, the winter precipitation is very low and the 
spring is windy and arid. Snow can only have very slight influence on peak biomass 
of the next year in Inner Mongolia grassland (Peng et al. 2010). Thus, the positive 
effect of snow is very weak in Ujimuqin grassland. However, we should point that 
the result of “no interaction effects” is just a region-specific finding for the Ujimuqin 
grassland, which is a wet part of the Mongolia plateau. Interaction effects of climate 
and density dependence may exist in the other regions of Mongolia plateau, especially 
in dryer grasslands. 

It should explain that why the population were density dependent in the three 
sites, while were not the others population. There is the large area of sandy land in 
sitel. The sandy land plant community is composed by shrub and tall grass with 
high biomass (Li et al. 2012). The tall plants are hardly completely buried by snow, 


and can provide some food during duzd. Thus, immigrating into sandy land is the 
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traditional adaption strategy of local herdsmen to against duzd. On another hand, 
there are many hills in the site6 and site7. The snow accumulation may be more 
heterogeneity; the snow depth may be deep in somewhere, and shallow in other place. 
We guessed that the landform and vegetation heterogeneity may be the reason why the 


population is density dependence in the three sites. 


The influence of climate and density dependence on nomadic society 

The impacts of climate and density dependence on domestic animal marked the 
behavior of nomadic herdsmen deeply. Although the nomadic Mongolian had a lot 
traditional ecological knowledge about duzd, but the only adaption strategy was 
migration (Jia 2011), which is very similar with the wild animal behavior (Stien et al. 
2010). The migration during pasture locked condition is full of risks for food lack 
and energy deficit (Weladji and Holand 2003). It recorded that sheep population can 
loss 30% individual during the five days migration in a duzd, this loss proportion of 
cattle and horse were about 20% (Wang 2006). During normal years, Mongolian 
herdsmen know it should reduce high animal density to avoid the outbreak of 
epidemic disease, thus they usually migrate out to keep low animal density. Many 
interviews of nomadic Mongolian people suggested that the migration distance and 
frequency positively correlate with the population size in nomadic period 
(Erdenebaatar 2003, Junichi 2005, Wang 2006, Jiang 2007). Therefore, just like 
these climate disasters, the density dependence effects can “push” the nomadic people 


to migrate out too. 
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Toynbee (1934) summarized that all nomadic people invasion are due to two 
external driving forces, one is the climate disasters in the steppe which push the 
nomadic people off the habitat, and the other is the social vacuum of neighborhood 
sedentary regions which pulls the nomadic people out the steppes. The “pushing 
force” and “pulling force” theory is succeed by most subsequent empirical studies 
after Huntington (1907) and Toynbee (1934); especially the “pushing of climate 
deterioration” are the core theory of explaining nomadic nation rise, invasion and 
conquest from natural environment perspective, and have derived many specific 
hypotheses, such as climate cooling driving hypothesis (Jenkins 1974, Lamb 1997), or 
cooling and drought driving hypothesis (Brown 2001). However, our findings 
suggest another “natural pushing force”-- density dependence, which would drive 
nomadic people to migrate when the animals are flourishing during climatic optimum. 
Therefore impacts of climate change on rise of nomadic nations may be more 
complex than precious understanding. Both bad climate and good climate can cause 
nomadic people migration. 

Especially, the “density dependence pushing force” theory can help us to fix the 
miss link of the “climatic optimum hypothesis” about the rise and conquest of 
Genghis Khan. It had already been known that the climate was warm and wet in 
Eurasian steppe during the Medieval Climate Optimum, and the harsh winter was less 
than other period of record history (HVIS 2011). Thus, the animal population would 
increase to a very high level, and then the density dependent effect would strengthen 


too. Under this scenario, nomadic people would migrate out to weaken the density 
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dependent effect, and then expand their habit. “The migration which was aim to 
weaken density dependence is more powerful than the migration for avoiding climatic 
disaster. During climate deterioration periods, the loss of animal and people would 
undermine the nomadic nation, and then it can not support an extensive conquest. 
Many nomadic empires resigned or collapsed after serious climatic disasters in history, 
such as Eastern Turk empire collapsed for a climate suddenly cooling in the 7th 
century (Fei et al. 2007), and Huns resigned to Han dynasty for drought. However, 
the density dependence can push nomadic people off the steppes with more human 

and animal recourse during climatic optimum. Now, the climatic optimum 
hypothesis seemed more credible to explain the Genghis Khan rise and conquest after 


it was assisted by “density dependence pushing force” theory. 


The impacts of natural regulation weakening in sedentary period 

Human behavior can not eliminate the natural factors influence on population 
dynamics in nomadic period, however neither climate nor population density can limit 
the population growth in sedentary period. The impacts of winter climate on animal 
are closed out the door, and the density dependent effects based on food competition 
and epidemic were eliminated by storing forage and modern veterinary knowledge. 
The simulation results showed the sheep population would still fluctuate and crash in 
the period of 1987-2005 if the negative feedbacks had not weakened. These 
simulation results indicated that the reason of overgrazing. 


Our study results may provide a new insight for future rangeland ecosystem 
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management in Inner Mongolia plateau. During nomadic period, both climate and 
density dependence are negative feedback which maintain the balance of animal 
population and vegetation dynamics. The rangeland ecosystems are spontaneously 
sustainable in this nomadic period. However, these natural feedbacks had been 
weakened or disappeared during sedentary period; and harvesting effect was the only 
limitation factor of animal population growth before the vegetation system collapse. 
The maintain mechanism in sedentary period was simpler than nomadic period, and 
more depended on human interventions (Fig. 7). Another study had found the 
vegetation dynamics were seldom influenced by grazing when the animal population 
was in the boom and bust cycle, while showed significant degeneration trend for the 
grazing pressure continue increased (Li et al. 2012). “Therefore, the rangeland 
ecosystem had shift from the natural spontaneity sustainable status into the human 
intervention dependent sustainable status (Fig. 7). Inner Mongolia people had 
gradually conquered the duzd, and achieved the flourishing of grazing husbandry, but 
the rangeland ecosystems resilience had been weakened meantime. This case 
suggested that the improvement of human adaption ability to natural disaster may 
cause unexpected ecosystem regime shift and result in more complex situation. Our 
study did not aim to blame the social change and effort of local herdsmen in 


combating with duzd, but we just try to understand how complex are the ecosystems. 


Conclusion 


In nomadic period, sheep population dynamic was regulated by combined effect of 
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climate and density dependence in Ujimuqin grassland of Inner Mongolia Plateau. 
There were several density dependence populations in the study area. For these 
density independence population, population growth mainly limited by population 
density in mild winters, but only controlled by climate during harsh winter years. 
These natural population regulation feedbacks limited the sheep population growth, 
and kept the population following the boom and bust cycle pattern. Because the 
density dependent effect is usually neglected in Mongolia plateau, our findings might 
be a “black swan” discovery. 

Our findings proved that Mongolian nomadic people would migrate for avoiding 
climatic disaster, or for weakening density dependence during the climatic optimum. 
Thus, we proposed there is “density dependence pushing force” for nomadic nation 
rise, invasion and conquest during climatic optimum. The “density dependence 
pushing force” changed the pervious common sense that only “climatic disaster 
pulling force” can cause nomadic population immigration. Besides, the “density 
dependence pulling force” view can fix the critical link of climatic optimum 
hypothesis about rise and conquest of Genghis Khan, and made it more credible than 
before. This finding of density dependence pulling force would help us to study 
more nomadic history in Eurasian steppe. 

On the other hand, the climate effect on population dynamics weakened, and the 
effect of negative density dependence disappeared during sedentary period. These 
population regulation mechanism changes had driven the sheep population continuing 


growth in the sedentary period. It indicated that the human society change can 


27 


weaken the natural population regulation mechanisms, and drive the rangeland 
ecosystems regime shift from “animal population fluctuation but ecosystem is 
spontaneously sustainable” status into “animal rapid increase but vegetation continues 
degeneration” status. “This knowledge of regime shift can provided a new 
explanation for recent overgrazing problem and rangeland ecosystem denegation in 


Inner Mongolia grassland. 
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Table 1. Bulmer R(R*) test results of the ten sites during the nomadic periods 


(1962-1979). 


If the statistics are smaller than critical value (a=0.05), then the Ho 


of is density independence rejected, and the population was marked as boldface. 


Statistics a=0.05 Site 1 Site 2 Site 3 Site 4 Site 5 
Bulmer R 0.8356 0.4540 3.9986 3.6842 0.9703 1.0578 
Bulmer R* - 0.4372 -0.1069 0.0439 0.0797 0.0404 0.1290 
Statistics a=0.05 Site 6 Site 7 Site 8 Site 9 Site 10 
Bulmer R 0.8356 0.9002 0.4703 0.9578 2.3731 1.4815 
Bulmer R* - 0.4372 -0.0400 -0.1637 0.0764 0.1220 0.0554 
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Table 2. The results of Pollard et. al’ test about 20 sheep population time series. 
The Ho of density independence was rejected when the 723 of original data is smaller 
than critical value of the 723* (a=0.05) which was generated from the random shuffled 


dataset. The density independence population was marked as boldface. 


Nomadic period (1962-1979) Sedentary period (1987-2005) 


Site T25*(a0.05) T23 of original T>3*(a=0.05) T23 of original 
data data 
1 0.7650 0.4607 0.7804 0.9623 
2 0.7703 0.9026 0.7803 0.8828 
3 0.7803 0.9484 0.7865 0.9994 
4 0.7672 0.7465 0.6895 0.6432 
5 0.7500 0.9357 0.7815 0.6330 
6 0.7664 0.6326 0.7880 0.9394 
7 0.7677 0.6150 0.7581 0.9999 
8 0.7695 0.7883 0.7790 0.9216 
9 0.7317 0.9880 0.7824 0.9179 
10 0.7681 0.7912 0.6906 0.9773 
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Table 3. The results of Dannis and Taper "test on 20 sheep population time series. 
The Ho of density independence was rejected when the T23? of original data is larger 
than critical value of the 725? * (a=0.05) which was generated from the bootstrap 


dataset. The density independence population was marked as boldface too. 


Nomadic period (1962-1979) Sedentary period (1987-2005) 


Sites T» (a0.05) T23 of Original T3 *(a-0.05) T23 of Original 
data data 
1 4.7840 16.5900 4.7040 0.5179 
2 4.9566 1.1407 4.4657 1.4947 
3 3.7077 0.5389 4.5440 0.0307 
4 4.8637 4.8536 5.1350 6.6623 
5 5.4164 0.7908 4.5338 9.8753 
6 4.9143 8.2752 4.6414 0.7821 
7 4.6947 10.4106 4.7122 0.0049 
8 4.8972 3.7811 4.4222 1.0985 
9 6.6599 0.0399 4.5227 1.2027 
10 5.5849 1.2576 4.7094 0.3054 


35 


Table 4. The optimal regular model and threshold model during nomadic period and 
was labeled as boldface. "v. “di”, “wi” and “K” stand for intrinsic growth rate, 
population density and winter index, parameter number of equation respectively. 


Site Method Optimal Models (1962-1979) K Adj-R? AICc 
Regular y, 21.28 — 0.31d, 2 0.51  -42.0 
1 
_ f1.37-0.34¢, if w, «1.69; 
Threshold yon 0.23-0.003w, jf w, 21.69. 4 0.51 -36.9 
Regular y; =0.35 —0.02d, “w, 2 016 -27.9 
2 
_ (0.83 - 0.154, if w, «171; a 
Threshold Ji 73975. 0.25w, f w217. 0.62 -29.0 
Regular y; = 0.42 - 0.18w; 2 0.33 -25.0 
3 
_ fo.18-0.144, if w, «142; 
Threshold 9i 730,84 —0.40w; P Wete. 4 0.77 -28.3 
Regular y; = 0.91 — 0.20d, —0.06w, 3 048  -340 
4 
_ [0.87 - 0214, if w, «1.94; 
Threshold Ji 7306 0.17w; f w2194. 4 0.61  -33.5 
Regular y; = 0.39 —0.03d, “w, 2 018 -24.9 
5 
_ [0.96 - 0.234, if w, «1.94; 
Threshold y= 0.90 — 0.31w; # w, 21.94. 4 0.62 -25.8 
Regular y; =0.94-0.24d; —0.06w, 3 0.57  -39.6 
6 
` f0.88-0.23d, if w «178; 
Threshold y= 119-031n, C w2178. 4 0.55 -36.5 
Regular y; =1.06-0.22d4; —0.04w, 3 0.42 -34.6 
7 
— [1.05 — 0.234, if w, «1.82; E 
Threshold 1 =70.63—0.1 pur if w218. 0.63 — -35.4 
Regular y; = 0.39 —0.03d, “w, 2 0.40 -32.1 
8 
_ f0.80-0.18d, if w, «115; 
Threshold JE e lan 4 0.77 -34.4 
Regular y; =0.43-0.21w, 2 034 -22.1 
9 
_ [0.50 - 0.084, if w, «175 
Threshold NIT 4 — 0,52w, # wèl. 4 0.68 -22.6 
Regular y; = 0.42 - 0.34d, “w, 2 0.39 -26.6 
10 
_ [1.07 — 0204, if w, <1.45; 
Threshold i. M scs 4 0.80 — -30.0 


36 


Table 5. The optimal models for each population in sedentary period (after 1987). 
The “yi”, “di”, “wi” and “K” are the same with table 4. 


Model 


Site Optimal Models K Adj-R? AICc 
Type 
3 y; =0.47—-0.07w, 2 0.25 -30.9 
8 y; = 0.45-0.04w, 2 0.29 -36.2 
Typel 
9 y; 2 0.49 — 0.06w, 2 0.14  -30.0 
10 y; 2 0.29 — 0.02w, 2 0.18  -494 
1 y, =-0.41+0.15d, 2 0.08 -34.1 
Type2 6 y; 2 0.17 - 0.10d, 2 0.12 -34.8 
4 y; =—0.09 + 0.17d, 2 0.51 -46.5 
5 y; 2 —0.14 4 0.19d, —0.04w, 3 0.53  -353 
Type3 2 y; =-0.20+0.18d; — 0.04w, 3 0.88 -31.6 


7 y; =-0.50+0.26d, — 0.03w, 3 049 -37.7 
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Figure Captions: 

Fig. 1. The geographic location of study area in the Eurasian steppe and population 
time series of this study. 

Fig. 2. “The demographic parameter (mortality, birthrate and harvested rate) from 
1962 to 2005 for each site. 

Fig. 3. The relationship between mortality, birthrate and winter climate during mild 
and harsh winters respectively at the regional scale. 

Fig.4 The relationship between mortality, birthrate and population density during 
mild and harsh winters respectively at the regional scale. 

Fig. 5. The regional relationship between intrinsic growth rate, population density 
and winter climate during mild or harsh winters respectively. The figure did not 
include the data of population at sitel, 6 and 7, for these population are density 
dependence. 

Fig. 6. The simulation population pattern of 1987-2005 under the nomadic scenario 
by the real climate condition of 1987-2005 and population abundant of 1987 at each 
site. 


Fig. 7. The rangeland ecosystem regime shift map. 
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Fig. 2 
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